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Abstract 

Dependence of magnetic field generation on the rotation rate is explored by direct numerical simulation 
of magnetohydrodynamic convective attractors in a plane layer of conducting fluid with square periodicity 
cells for the Taylor number varied from zero to 2000, for which the convective fluid motion halts (other 
parameters of the system are fixed). We observe 5 types of hydrodynamic (amagnetic) attractors: two 
families of two-dimensional (i.e. depending on two spatial variables) rolls parallel to sides of periodicity 
boxes of different widths and parallel to the diagonal, travelling waves and three-dimensional "wavy" rolls. 
All types of attractors, except for one family of rolls, are capable of kinematic magnetic field generation. 
We have found 21 distinct nonlinear convective MHD attractors (13 steady states and 8 periodic regimes) 
and identified bifurcations in which they emerge. In addition, we have observed a family of periodic, two- 
frequency quasiperiodic and chaotic regimes, as well as an incomplete Feigenbaum period doubling sequence 
of bifurcations of a torus followed by a chaotic regime and subsequently by a torus with 1/3 of the cascade 
frequency. The system is highly symmetric. We have found two novel global bifurcations reminiscent of the 
SNIC bifurcation, which are only possible in the presence of symmetries. The universally accepted paradigm, 
whereby an increase of the rotation rate below a certain level is beneficial for magnetic field generation, while 
a further increase inhibits it (and halts the motion of fiuid on continuing the increase) remains unaltered, 
but we demonstrate that this "large-scale" picture lacks many significant details. 

Keywords: Rayleigh-Benard convection, convection in rotating fluid, kinematic dynamo, nonlinear 
magnetohydrodynamic regimes, bifurcations 
PACS: 47.20.Bp, 47.20.Ky, 91.25.Cw 



1. Introduction 

Magnetic held of stars, planets and other astrophysical objects is often attributed to the motion of 
electrically conducting melted substance in their interior (Parker, 1979; Priest, 1984; Soward et ai, 2005; 
Hughes et ai, 2007; Dormy and Soward 2007), which is usually sustained by compositional and thermal 
convection. Convective flows in a plane layer, ranging from very simple (Matthews, 1999) to turbulent 
(Meneguzzi and Pouquet, 1989; Cattaneo et ai 2003) ones, are capable of magnetic field generation. In 
simulations of the dynamo in the Earth's liquid core (Glatzmaier and Roberts, 1995; Roberts and Glatzmaier, 
2000), a magnetic field of the approximately correct strength was produced, which had the dipole structure 
and exhibited reversals similar to the natural ones. Dynamos in spherical shells were also simulated by 
Grote and Busse (2001), Ishinara and Kida (2002), Takahashi and Matsushima (2005) and other authors. 
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Magnetohydrodynamic thermal convection is characterised in the dimensionless form by the Rayleigh 

number, R (indicating the magnitude of thermal buoyancy forces); the Prandtl number, P (the ratio of 
kinematic viscosity to thermal diffusivity) , the magnetic Prandtl number, (the ratio of kinematic viscosity 
to magnetic diffusivity), and the Taylor number, Ta (measuring the speed of rotation). Usually no-slip or 
stress-free boundary conditions for the flow velocity and perfectly conducting or insulating boundaries for 
magnetic field are considered. 

How magnetic field generation by convection depends on parameter values, is not yet explored in full. 
Large Rayleigh numbers arc beneficial for generation (the critical Pm decreases monotonically with the 
increasing R) in spherical shells (Bussc, 2000). Podvigina (2006) found a similar dependence in plane layer 
dynamos on increasing R over the critical value, but for R over a certain threshold the behaviour of the 
critical Pm ceased to be monotonic; Podvigina (2008) also examined the influence of the Prandtl number. 

Rotation, which is a common feature of the majority of astrophysical bodies, can also assist magnetic 
field generation. Rotation of the Earth is relatively rapid; it is believed, that geodynamo operates in the 
outer core in the magnetostrophic regime, in which the strength of the Coriolis force is comparable to that 
of other primary forces - the Lorentz force, pressure and buoyancy. Rotation is also an important factor in 
physics of the Solar tachoclinc (Christensen-Dalsgaard and Thompson, 2007), where apparently the solar 
dynamo is located. Differential rotation gives rise to cj— effect dynamos (Moffatt, 1978). Although they are 
slow, as opposed to fast ones, which are believed to operate in stars, this mechanism is regarded as a key 
element of the dynamo of the Sun (Tobias and Weiss, 2007). Boundary layers and shear flows developing in 
rapidly rotating fluids are the structures controlling the dynamics of fluid and planetary dynamo processes 
(Busse et ai, 2007). Therefore, how the processes of generation are affected by the rate of rotation is an 
interesting question for astrophysical applications. 

Dependence of magnetic field generation on the rate of rotation was explored in a number of papers. 
Meneguzzi and Pouquet (1989) observed that rotation could be beneficial for nonlinear dynamos - the 
critical magnetic Reynolds number decreased and the ratio of magnetic to kinetic energies significantly 
increased, when rotation was on. By contrast, Cattaneo and Hughes (2006) found that rotation was not a 
significant factor: they reported "similar growth rates and similar saturation levels" in rotating and non- 
rotating systems. Only several runs were presented in eac;h of the two papers, and hence the observations 
were inconclusive. Whilst turbulent convective nonlinear dynamos do not require rotation, near the onset of 
convection in a plane layer rotation is essential for both kinematic (Matthews, 1999) and nonlinear dynamo 
action (Demircan and Seehafer, 2002) - in both regimes considered by these authors dynamos failed in the 
absence of rotation. 

Our paper is devoted to investigation of the dependence of magnetic field generation on To. A fluid 
heated from below in a plane horizontal layer rotating about the vertical axis is considered in the Boussinesq 
approximation, whereby the buoyancy depends linearly on temperature, density variation is neglected in 
the mass conservation equation and the flow is incompressible. Perfectly electrically conducting stress- free 
horizontal boundaries of the layer are held at constant temperatures; periodicity in horizontal directions 
with the same period L (measured in the units of the layer depth) is assumed. 

We flx all parameter values except for the Taylor number and investigate numerically attractors of 
the system for Ta increasing from zero, and bifurcations delimiting branches of the attractors. Along 
the branches we trace average magnetic, E^, and kinetic, E/-, energies, as well as their ratio, Em/Eh, in 
saturated regimes; the latter quantity is a measure of dynamo efficiency of prime concern in astrophysics 
(although we arc clearly not in the astrophysical range of parameter values). The system is also of interest 
from the point of view of equivariant bifurcation theory, because it has a large symmetry group. 

Which parameter values are optimal for such an investigation? Naturally, we need to employ the values, 
for which the structure of attractors of the dynamical system and the geometry of their branches in the 
parameter space are relatively simple. Behaviour of the dynamical system, which one encounters in the geo- 
and astrophysical environments or in experimental dynamos, results from a large number of bifurcations of 
the trivial steady state, and it is prohibitively complex to serve as a starting point for such a study. These 
difficult cases can be approached later by continuation in parameters starting at the regimes determined in 
the initial study. Also, in the initial investigation we cannot employ relatively small values of Pm, because 
for them only vigorous turbulent flows would yield magnetic dynamos: Many hydrodynamic bifurcations 
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would be required to bring the convective MHD system into the turbulent state and, consequently, for such 

parameter values the system will also be prohibitively complex. By contrast, for relatively large Pm already 
laminar convective flows of simple structure can act as dynamos; moreover, if Pm is not too large, magnetic 
field will not be too strong to enable the Lorentz force to destabilise the convective flows. This suits ideally 
our purposes. 

Following these considerations, we employed the same values as in Podvigina (2006): P = 1, Pm = 8, 
R = 2300 and L = 2^/2. The Rayleigh number is not far from the critical value for the onset of convection; 
hence, convective attractors have a simple roll structure, which simplifies a detailed investigation of the 
attractors and bifurcations. The Pm value is not far from its critical values for the kinematic dynamo problem 
for convective attractors, hence flows in convective MHD attractors (when dynamos operate) are qualitatively 
similar to flows in non-magnetic ones. (Note, that in steady MHD states flows can be reconstructed from the 
structure of the generated magnetic field, see Zheligovsky, 2009a.) The assumed aspect ratio, L, is equal to 
the spatial period (measured in the units of the layer width) of the hydrodynamic mode becoming unstable 
the first when the onset of convective motion occurs in a non-rotating layer on increasing the Rayleigh 
number. 

Although these values are chosen on the basis of mathematical convenience, they arc not unphysical: 
Pm are of the order of unity in accretion disks (von Rckowski et al., 2001): magnetic Prandtl numbers 
Pm in excess of unity arc typical for intergalactic and interstellar gases (Shukurov and Sokoloff, 2008; see 
also Table 1 in Brandenburg and Subramanian, 2005). A plasma experiment for the study of astrophysical 
dynamos, where 3-10^^ < P„i < 5, was proposed by Spence et a/.(2009). We also note that typically galaxies 
are thin disks (see Ruzmaikin et al., 1988), and hence some features of the problem at hand (an infinite 
plane layer with stress-free boundaries) are appropriate to model galactic dynamos. We stress nevertheless, 
that we solve the problem in an abstract setting and do not pretend to simulate any specific physical system. 



2. Statement of the problem 

The system is governed by the Navier-Stokes equation 

— = V X (V X v) -h PAv + PRee^ + Ptv x e^ - Vp - b x (V x b), (l.a) 
the magnetic induction equation 

^ = Vx(vxb) + PP-iAb (1.6) 

for solenoidal fields 

V-v = 0, V-b = 0, (l.c) 

and the heat transfer equation 

— = -{,r-S7)0 + v,+Ae. (l.d) 

Here v denotes the flow velocity, b the magnetic fleld, 6 the difference between the temperature of fluid and 
the linear temperature profile, and r = y/Ta is twice the angular speed of rotation of the fluid. 
The following boundary conditions on the horizontal boundaries: 

^ = ^ = v^=0, 9 = at^ = 0,l (2.a) 
oz oz 

^ = ^ = 6. = at^ = 0,l (2.6) 
and periodicity in horizontal directions with the same period 

v{x,y, z) ~ v{x + mL,y + nL,z), 6{x,y, z) ^ 6{x + mL,y + nL,z), (3. a) 

h{x,y, z) =h{x + rnLjy + nL, z) (3.6) 
Vm, n e Z 

are assumed. 
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The equations are solved numerically by the standard pseudospectral methods (Boyd, 2001; Peyret, 
2002). Fields are represented as Fourier series satisfying the boundary conditions (2): 

/ -0^cos(7rn3z) \ 
v = ^ *nCOs(7rn3z) e^("i^+"^^), 
n y '0^sin(7m3z) / 

(6^cos(7rn3z) \ 
6^cos(7rn3z) e'^'("i"+"^^), 
6^sin(7rn3z) / 

e^J2^n sin(7rn3z)e^("i^+"^^). 

n 

The resolution of 31 x 31 x 17 Fourier harmonics has been employed for computation of hydrodynamic 
convective attractors, and 63 x 63 x 33 in simulations of kinematic dynamos (without dealiasing in both 
cases). These simulations have been checked against runs with the double resolution without dealiasing, and 
computations of kinematic magnetic modes also against runs with the resolution of 41 x 41 x 21 harmonics 
with dealiasing. We have been employing the standard method of dealiasing in computation of all products 
in (1) (requiring to evaluate fields on a uniform 64 x 64 x 33 mesh in the physical space in order to compute 
41 X 41 X 21 Fourier harmonics of the products). 

Simulation of nonlinear hydromagnetic convective regimes has been performed with the resolution of 
41 X 41 X 21 Fourier harmonics with dealiasing. We note that computation with a coarser resolution of 
31 X 31 X 16 Fourier harmonics without dealiasing in earlier MHD convection simulations by Gertsenshtein 
et al. (2007, 2008) yielded a wrong classification of some regimes (a typical error consisted of obtaining a 
periodic regime instead of a convective MHD steady state). Although from the conservative point of view 
it thus may be also unsafe to employ the resolution of41x41x21 harmonics, the huge amount of runs 
has forced us to use it in the main bulk of computations. For this resolution magnetic energy spectrum 
decreasef0 by at least four orders of magnitude and the flow and temperature energy spectra decrease by 
about nine orders of magnitude. Some computations of nonlinear hydromagnetic regimes (including at least 
one attractor on each MHD branch) have been checked against runs with the resolution of 127 x 127 x 65 
Fourier harmonics without dealiasing, the results remaining visibly unaffected (the energies being reproduced 
with the accuracy better than 0.001%). 

By Ek and Em we denote kinetic and magnetic energies, respectively, averaged over the fluid layer: 

1 IvP 1 IbP 

^k^jlj J J ^rfa^rfyc?^, Ern^j^J J J ^-^dxdydz. 

For non-steady convective MHD regimes, time averaging is also performed. 

Branches of attractors are traced by continuation in parameter: computations are done for initial con- 
ditions, which are a point (in the phase space) on the attractor for a "nearby" Ta. Typically the distances 
between such "nearby" Ta vary between 0.1 and 100. Also, for some Ta runs have been performed for 
"random" initial conditions, comprised of fields with pseudorandomly generated Fourier coefficients and an 
exponentially decaying spectrum, with either small (^ 10~^), or large (~ 10^ — 10^) initial kinetic, magnetic 
and thermal energies, in order to check whether multiple attractors coexist for the value of the Taylor num- 
ber under consideration. Bifurcations of steady states are located by solving the eigenvalue problem near 
the endpoints of the branches and extrapolating the eigenvalues or their real parts to zero. 



^The ratio of energy in the spherical shell of width one in the Fourier space with the largest energy content to the energy 
in the last considered spherical shell is reported. 
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3. Symmetries 

The symmetries of the rotating hydromagnetic convective system under consideration constitute a sub- 
group of the group of symmetries of the system in the absence of rotation. We hst them here for reader's 
convenience, using the notation of Podvigina (2006). 

The symmetry group of the convective system (1) with the boundary conditions (2), (3) is Z4 K x Z2. 
The group Z4 consists of rotations 

si : {x,y,z) H> {y,-x,z), 

52 ■■ ix,y,z) {-x,-y,z), 

53 : {x,y,z) H> {-y,x,z) 

and the identity sq = e. T^^ and Ty are the groups of translations in the x and y directions, respectively: 

la ■ ix,y,z) i~> {x + a,y,z), 
la ■ ix,y,z) ^ {x,y + a,z) 

where < a < L (7£ =7^ = c)- "^xy is the group of translations along the diagonal: 

la^ ■ {x,y,z)>-^ {x + a,y + a,z). 

The group Z2 is generated by reflections about the horizontal midplane: 

r : {x,y,z) 1-^ {x,y, 1 - z). 

If magnetic field is present, the group of symmetries of the system is augmented by the symmetry 
reversing magnetic fleld 

q: (v,^,b)^(v,^,-b). 

The symmetry about a vertical axis, S2, and parity invariance about a centre (located on the midplane) 
which is a composition of S2 and r, are of special interest, since in such MHD states the global a— effect is 
zero: In the presence of the global a— effect the system is inherently unstable to large-scale perturbations, 
while in its absence the instability, if present, develops on time scales of a higher order (see Zlieligovsky, 
2009b). Although we have not studied spatio-temporal symmetries of non-steady attractors in detail, we note 
that we have encountered periodic orbits possessing the symmetry about a vertical axis or parity invariance 
with a time shift (equal to a half of the temporal period) . A spatio-temporal symmetry with a time shift 
differs from the respective spatial symmetry in that it relates fields at two time instances, separated by an 
interval whose length is equal to the time shift; for instance, the symmetry S2 about a vertical axis through 
the point (a^, ay,0) with a time shift T of a vector fleld f is deflned by the conditions 

fx{ax -x,ay -y,z,t) = -fx{ax+x,ay+y,z,t + T), 

fy{ax -x,ay - y,z,t) = -fyia^ + x,ay + y,z,t + T), 
fzidx -x,ay- y,z,t) = f^{ax + x,ay + y,z,t + T). 

Note that compositions S27a ^^'^ *2''7a are the symmetry about a vertical axis and parity invariance, in 
which the axis (the centre, respectively) of symmetry is displaced in the direction d by —a/ 2; changing the 
order of factors in these compositions to the opposite one results in the displacement of the axis and the 
centre, respectively, in the reverse direction. S2 is used as a generic notation for the symmetry about a 
vertical axis with no specific axis indicated; consequently, we do not distinguish notationally §27^ from S2 
(except for in Table 5). In what follows, for construction of figures displaying vector fields which possess 
the symmetry S2 we employ a coordinate system with the origin lying on the axis of 82- 
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Figure 1: Kinetic energy (vertical axis) of liydrodynamic attractors for < Ta < 2000 (horizontal axis). Labelling of attractors 
is explained in Section 4 (see also Table Note that Ri does not bifurcate from Rd: energies of the two flows at the point, 
where Ri becomes unstable, differ by about 0.5, which is not seen in the scale of the figure. 

4. Hydrodynamic convective attractors 

For a small Rayleigh number, i.e. a small temperature difFerence between the upper and lower boundaries, 
the fluid is not moving and heat is transported by thermal diffusion only. When R exceeds the critical value 
(which increases with Ta), a fluid motion sets in. For P = 1 considered here, the instability is monotonic, the 
primary motion has the form of rolls. The wave number (in the horizontal direction) of the most unstable 
mode monotonically increases with Ta. If periodicity in horizontal directions is imposed, the horizontal 
wave number must be compatible with the periodicity cell size and it can take only discrete values, and for 
an increasing Ta the wave numbers of the emerging rolls constitute a piecewise constant growing function. 
At very high rotation rates an array of very thin rolls emerges (Bassom and Zhang, 1994). 

Five types of attractors of hydrodynamic convection (governed by the system (l)-(3) for b = 0), that 
we have found in computations, are shown on the bifurcation diagram Fig. [T] (see also Table [T|). Since the 
Rayleigh number is not far from the critical value for the onset of convection, the flows are of a simple 
spatial structure (see Fig. 

Three types of attractors are steady rolls of different spatial periods. In agreement with the linear 
stability theory, their wave numbers increase with Ta (remaining the same within each branch). The rolls 
obtained in computations are parallel to a side of the periodicity cell with either the same period as the 
periodicity cell size, or a half of it; or they are parallel to the diagonal of the cell. The branches are labelled 
Ri, R2 and Rd, respectively. 

In a non-rotating layer, bifurcations of rolls are well documented. However, these results cannot be 
directly applied to convection in a rotating layer, even if the rotation rate is small - this can only be done 
after unfolding of bifurcations of rolls is performed, i.e. after it is determined how the bifurcations alter, if 
small terms breaking the reflection symmetry are added. This is an interesting problem in its own right, 
but it is beyond the scope of the paper. 

On increasing R, in a non-rotating layer rolls can bifurcate to travelling waves in a Hopf bifurcation 
for P < 1 (Getling, 1998). This bifurcation, with the 0(2) symmetry group, was studied in detail in 
Golubitsky et al. (1988). In it, branches of standing and travelling waves emerge; if both branches bifurcate 
supercritically, one of them is stable. A travelling wave, TW, is time-periodic in a coordinate frame at 
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Table 1: Attractors found numerically for hydrodynamic (b = 0) convection. Column 2 presents the interval, where existence 
of branches of attractors is confirmed numerically, column 3 presents the symmetry group for which an attractor is pointwise 
invariant, column 4 generators of the group (for an appropriately chosen Cartesian coordinate system); if a group is a product 
of several subgroups, generators of the subgroups are separated by semicolons. Column 5 presents locations (Ta) and types of 
bifurcations in which a steady state becomes unstable (S denotes a steady-state bifurcation, H a Hopf bifurcation), column 6 
the dimension of the respective center eigenspace, column 7 the action of the symmetry group on the null eigenspace, and the 
last column elements of the group which act trivially. 1 denotes a trivial group of symmetries. 



Type 


Interval of 


Symmetry 


Generators 


Bifurcation 


D 


Action 


Kernel 




existence (Ta) 


group 












TW 


[1,187] 


Z2 












Ri 


[188, 725] 


D2 K T 




187.19, H 


4 


0(2) 


nl/2 








725.45, S 


2 


0(2) 


nl/2 


WR 


[719, 750] 


D2 




718.16, S 


1 


1 


nl/2^ S2 










750.82, S 


1 


Z2 


s^nl/2 


Rd 


[684, 1373] 


D2 X T 




683.64, S 


2 


0(2) 


nl,2 










1373.32, S 


2 


0(2) 


x,—y 

nL/i 


R2 


[999, 1969] 


D4 K T 


nl/i. S2; 7^ 


998.18, S 


2 


0(2) 


xy 

nL/4 



rest and it is steady in a frame moving with the speed of the pattern. In a non-rotating convection such 
bifurcation from Ri to a stable travelling wave takes place near R « 1755 (Podvigina, 2006). The instability 
causes a sinusoidal bending of rolls, the pattern travelling along the axis of a roll. 

We also observe the TW: On decreasing Ta, Ri loses stability and the stable travelling wave emerges 
in a similar supercritical Hopf bifurcation with the 0(2) symmetry group. Due to rotation, no reflection 
symmetries are present, enabling the travelling wave to drift in the direction perpendicular to axes of rolls. 
All other bifurcations of rolls of the three types are subcritical pitchfork bifurcations, with no stable branches 
emerging. On increasing Ta, the branch R2 terminates at Ta = 1969.67 (this number is in a good agreement 
with the value of the Taylor number, obtained from the Chandrasekhar's (1961) formula for the critical value 
of the Rayleigh number for the onset of convection in a layer with free boundaries) on the trivial steady 
state; no non-trivial convective regimes exist for higher rotation rates. 

Another attractor, distinct from rolls, is wavy rolls, WR. In a non-rotating layer emergence of WR from 
Ri is a consequence of the 1 : ^/2 mode interaction (Podvigina and Ashwin, 2007). Similar arguments 
show that in a rotating layer WR can also bifurcate from rolls. On decreasing Ta, WR disappear in a 
saddle-node bifurcation - the branch turns back and becomes unstable. Apparently, it bifurcates from Ri 
at Ta = 725.45: WR and Ri have similar shapes and close energy values, and by the equivariant branching 
lemma (Golubitsky et at, 1988) the action of the symmetry group of Ri (see Table[T]) implies that a branch 
with the symmetries of WR bifurcates from Ri . 

5. Magnetic field generation 

We have investigated whether the flows that are convective hydrodynamic attractors are capable of 
kinematic magnetic field generation. For steady hydrodynamic states we have been computing (applying the 
algorithm of Zheligovsky, 1993a, b) dominant eigenvalues of the magnetic induction operator. For travelling 
waves, equation (2.b) in the co-moving reference frame yields a similar eigenvalue problem (see Podvigina, 
2006). 

All the hydrodynamic attractors have non-trivial symmetry groups and magnetic modes can be classified 
in the terms of the action of their symmetries. Symmetry groups of rolls are continuous, they include shifts 
along the axis of rolls, and bifurcating modes can have arbitrary periods in this direction. We restrict our 
attention to magnetic modes that have the periodicity of the hydrodynamic convective attractors. The 
computed growth rates are shown on Fig.[3l and symmetries of dominant modes are presented in Table [2l It 
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Figure 2: Isolines (step 2) of i;^ on the horizontal midplane z = 1/2 for hydrodynamic attractors TW, Ta = 50 (a); Ri, 
Ta = 500 (b); WR, Ta = 720 (c); Rd, Ta = 685 (d); R2, Ta = 1100 (e). Solid lines indicate non-negative, dashed lines 
negative values, x: horizontal axis, y: vertical axis. 
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Figure 3: Growth rates (vertical axis) of magnetic field generated kinematically by hydrodynamic attractors as a function of 
Ta (horizontal axis). 
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Table 2: Dominant kinematic magnetic modes of the hydrodynamic convective attractors. Column 2 presents generators for 
symmetry groups of the hydrodynamic attractors, column 3 the interval of Ta where a magnetic mode is dominant, column 
4 dimension of the eigenspace associated with the dominant eigenvalue, column 5 the action of the symmetry group on the 
eigenspace, column 6 the symmetries which act trivially, the last column the eigenvalue A of the magnetic induction operator 
for which the maximal growth rate is attained in the interval of Ta specified in column 3. 



Flow 


Generators 


Ta 


D 


Action 


Kernel 


X 


TW 




[1,84] 
[85, 187] 


2 
2 


Z2 
Z2 


inl/2 
nl/2 


0.97 ±28i 
1.18 ± 19i 


Ri 




[188,725] 


2 


0(2) 


nl/2^ 97l/2 


1.53 


WR 


S2, ni/^ 


[719, 750] 


1 


Z2 


9^2, »-7£/2 


0.45 


Rd 


S2, ''7E/2; 7"^^ 


[684, 925] 
[930, 1373] 


2 
2 


0(2) 
0(2) 


^97£/2. 97l/4 
^97l/2. Hl% 


0.49 
0.13 


R2 


S2, nl/i, 7^ 


[999, 1420] 
[1430, 1969] 


2 
2 


0(2) 
0(2) 


nl/4, qil/2 

xy y 


-0.48 
-0.62 



turns out that all the attractors, except for R2, can generate magnetic field. Note, that since TW bifurcates 
from Ri , the steady magnetic mode of Ri becomes a time-periodic mode of TW, whose frequency at the 
point of bifurcation coincides with the one of TW and the eigenvalue becomes the real part of the TW 
magnetic eigenvalue. 

There is no obvious relation between the structure of magnetic field generated by TW and stagnation 
(in the co-moving reference frame) points of the flow (Podvigina, 2006, found that the same was true for 
convective dynamos without rotation). The spatial structure of growing magnetic modes is shown on Fig. 2] 
and O For all convective flows capable of magnetic field generation, in dominant magnetic modes the field 
concentrates near horizontal boundaries in flattened half-ropes. A plausible underlying physical mechanism 
for this kind of behaviour in the case of perfectly conducting boundaries was proposed by St Pierre (1993). 
For steady flows, in agreement with the kinematic dynamo theory (Galloway and Zheligovsky, 1994), each 
half-rope is centered at a stagnation point of the flow and oriented along the one-dimensional unstable 
manifold of the flow at this point (see Fig. [5] (a), (c) and (e) ). Magnetic field is advected by the flow and, 
accordingly, the ropes are stretched along the trajectories of fluid particles on the upper boundary; however, 
advection is affected by magnetic diffusion, which must be responsible for a (rather modest) deviation of 
the direction of half- ropes near their ends from the direction of the trajectories (see Fig. |6] (a), (c) ). The 
half-ropes extend till they begin to feel the influence of the adjacent stagnation points in the direction of 
their stretching, which results in termination of the ropes. Magnetic fleld is advected by the flow into the 
layer; it is redistributed to form vertical two-dimensional magnetic flux sheets (shown on Fig. [S] (d) and (f ) ) 
in the plane of the unstable/neutral directions of the stagnation points. 

Sample trajectories of fluid particles inside the rolls, also shown on Fig. [6] (a), (c), (e), attest that 
the motion is not planar (and hence the Zeldovich, 1956, antidynamo theorem is unapplicable). In Ri 
and R2 the fluid moves along closed loops, because the flows are two-dimensional (i.e. independent of a 
horizontal Cartesian coordinate, which in our notation is y) and possess two symmetries, S2 (more precisely, 
the symmetry about any vertical line on the boundary between two adjacent rolls, rotating in opposite 
directions) and r^l/2- (Due to solenoidality and two-dimensionality, and Vz are associated with a stream 
function, say, iplxjz), and equations of motion for the particles imply that their trajectories are spirals 
residing on the surfaces tpix, z) = const, topologically equivalent to infinite cylinders. The two symmetries 
of the flow imply ip{x, z) = 'ipi—x, —z) and Vy{x, z) = ~Vy{~x, —z), in a coordinate system with the origin on 
the axis of the roll; hence the helical trajectories degenerate into closed loops.) The structure of trajectories 
in WR is by far more complex (two trajectories in the interior of the volume of fluid are shown on Fig. |5] 
(c) ), however, visibly this does not affect the complexity of the magnetic field ~ neither at large magnetic 
energy levels, nor at small ones. 
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(g) (h) 

Figure 5: Isosurfaces of magnetic energy density of the dominant magnetic modes, at the level of a half of the maximum, for 
Ta = 50, TW (a); Ta = 150, TW (b); Ta = 500, Ri (c); Ta = 720, WR (d); Ta = 685, Rd (e); Ta = 1200, Rd (f); Ta = 1100, 
R2 (g); Ta = 1700, R2 (h). Isolated stagnation points and lines of stagnation points of the flow on the horizontal boundaries 
are shown by dots and bold lines. Four periodicity cells are displayed. On each panel, coordinate axes are as shown on (a). 
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(a) (b) 




(e) (f) 

Figure 6: Isosurfaces of magnetic energy density of tlic dominant kinematic magnetic modes generated by Ri for Ta = 500 at 
the 50% level of the maximum (a) and at 3.5% (b); by WR for Ta = 720 at 50% (c) and 10% (d); by R2 at Ta = 1700 at 50% 
(e) and 12% (f). Stable and unstable directions of some stagnation points are shown by lines with arrows, neutral directions 
associated with zero and imaginary eigenvalues by lines without arrows. (See location of isolated stagnation points of the flow 
on the horizontal boundaries and lines of stagnation points on Fig. [5] and note that axes of rolls Ri and R2 are also comprised 
of stagnation points.) Sample trajectories of fluid particles on the upper boundary and in the interior of the fluid layer are 
shown by thin lines. One periodicity cell is displayed. On each panel, coordinate axes are as shown on Fig. [5] (a). 
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Fig. HI (b), (d), (f), showing isosurfaces of magnetic energy density at low levels, expose two-dimensional 
structures of magnetic field apparently related to stagnation points (for which eigenvectors of the stress 
tensor \\dvx,y,zld{x,y,z}\\ are also shown) of the flows. Formation of a magnetic flux sheet tangent to the 
two-dimensional unstable manifold of a stagnation point of a flow was investigated analytically by Childress 
(1979) and Childress and Soward (1985). The structures that we observe are peculiar in that magnetic 
energy density increases not in between the isosurfaces, as it does in magnetic flux sheets, but, on the 
contrary, outside them, i.e. the isosurfaces reveal two-dimensional "magnetic flux sheet gaps" . 

At first glance, emergence of such "antistructures" can be linked to the fact that in most cases the 
respective stagnation points have neutral eigendirections along lines of stagnation points (note that lines of 
stagnation points reside not only on the horizontal boundaries, as shown on Fig. [5l but also constitute the 
axes of the rolls Ri and R2) - at such stagnation points the two-dimensional antistructures are associated 
with the unstable and neutral eigendirection. In the presence of a neutral direction the theory of Childress 
(1979) and Childress and Soward (1985) may be unapplicable. This argument, however, breaks for isolated 
stagnation points of WR in the middle of upper edges parallel to the y axis (see Fig. El (d) ) , which do 
not possess neutral eigendirections, but have genuine two-dimensional unstable manifolds. By contrast, 
formation of "bells of trombones", the antistructures on Fig. [5] (b), is not related to any unstable direction 
- they are oriented along the eigenplane associated with two imaginary eigenvalues of the stress tensor. 
Dynamics of fluid in this case is significantly different from that considered ibid.: the trajectories are closed 
oval loops centered at the axis of the roll. Fig. H] (c) reveals, that the "bells of trombones" constitute 
boundaries between blobs of magnetic field of the opposite orientation. 

6. Nonlinear magnetic field generation 

Kinematic magnetic growth rates are not large (see Fig. In agreement with this, for all convective 
MHD attractors that we have computed, magnetic energy, Em, is smaller (at least seven times) than the 
kinetic one, Ek (see Figs. [THS]). Consequently, all the MHD attractors can be regarded as perturbations 
of hydrodynamic attractors (note that the plot of the ratio Em/ E^ is quite similar in shape to the plot of 
magnetic energy, cf. Figs. [3] and [51 also suggesting that we are close to the onset of both convective motions 
and of magnetic field generation) and the spatial structure of the flows is similar to that in the absence of 
magnetic field (cf. Figs. and [TU)) . Hence, we label branches of MHD attractors as follows: The main label 
denotes the type of an attractor: S for a steady state, P for periodic, Q for quasiperiodic (in assigning these 
labels we disregard temporal periodicities due to drifts along the x and y axes). The superscript denotes 
the hydrodynamic attractor, the MHD attractor is genetically related to. The subscript is the consecutive 
number within the collection of attractors of the specified morphology (e.g., there are eight different MHD 
steady states with spatial structure similar to that of Ri). The numbering of branches is in the order of 
increasing Ta. 

Despite hydrodynamic attractors of just four types generate magnetic field in the kinematic regime and 
magnetic field growth rates are not large, we have observed a large variety of MHD attractors (see Table [3] 
and bifurcation diagrams Figs. [7] and [5]). In the description of branches of convective MHD attractors and 
bifurcations bounding them we follow the ordering of hydrodynamic attractors in Table [TJ and for each of 
them we start with unstable magnetic modes. When on increasing P,n a magnetic mode eigenvalue crosses 
the imaginary axis, an MHD steady state or periodic orbit appears; for supercritical bifurcations these objects 
are stable. We refer to them as primary magnetic attractors. For Pm — 8, for which the problem is solved 
here, these objects are not necessarily stable; for a varying Ta they constitute a branch, and somewhere along 
the branch can gain stability. Because of emergence of the branches of unstable convective MHD states, 
in the present problem identification of hydrodynamic attractors with neutral kinematic magnetic modes, 
giving rise to the branches of primary convective MHD attractors, is not straightforward. We perform 
this identification by comparing the spatial structure of the fields in the attractor with the respective 
hydrodynamic attractor and the magnetic mode, including their symmetry groups. 

To identify bifurcations of steady states we calculate eigenvalues and the associated eigenspaces of the 
operator of linearisation of the system (l)-(3). Hence for each bifurcation of a steady state we know the 
dimension and the action of the steady state symmetry group on the eigenspace (see Table |4]) , which is a 



13 




14 




15 





(a) 

0.03 



800 



1000 1200 1400 



0.00 



\cRl 






,cRl V8 

Npe \ 












V \\pf 




gWR 





216 218 220 222 224 650 
(b) 



'00 750 800 

(c) 



Figure 9: The ratio of energies Em/E)^ (vertical axis) of MHD attractors for < Ta < 2000 (horizontal axis) 
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Table 3: MHD attractors found in computations. Column 3 presents the symmetry group for which an attra<;tor is pointwise 
invariant, column 4 generators of the group (for appropriately chosen location of the origin of the coordinate system); if a 
group is a product of several subgroups, generators of the subgroups are separated by semicolons. Columns 5 and 6 present 
time-averaged kinetic and magnetic energies, respectively. 



Type 
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existence (Ta) 
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Table 4: Bifurcations of MHD steady states. Column 2 presents the symmetry group of a steady state, column 3 the generators 
of the symmetry group, column 4 the critical Ta, column 5 the type of a bifurcation (S.P. denotes subcritical pitchfork, S.H. 
subcritical Hopf and S. saddle-node bifurcations), column 6 the dimension of the respective center eigenspa<;e, column 7 the 
action of the system symmetry group on the eigenspace, and the last column elements of the group which act trivially. Note 
that S2 mentioned in columns 3 and 8 have the same axis. 
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necessary information for application of the general theory of bifurcations for symmetric systems (Golubitsky 
et ai, 1988). 

For each growing magnetic mode there exists a primary MHD attractor. At points of bifurcations from 
rolls, dimension of the kernel of the magnetic induction operator is two due to the presence of translation 
symmetries (see Tabled]). The action of the symmetry group on the eigenspace is 0(2) generated by S2 and 
translations. Along axes of rolls eigenmodes can have periods £/n with an integer n > 0, where for rolls 
parallel to coordinate axes £ — L and for rolls parallel to a diagonal i = \plL. If a magnetic mode has a 
period £/n, the shift by £/{2n) acts as — /, hence the superposition of the shift by l/{2n) and q maps the 
mode into itself. Bifurcations from rolls are pitchfork with the symmetry group 0(2). Thus a continuum 
of steady states emerges, the symmetry group of each of them is a product of Z2 generated by S2 and the 
subgroup which acts trivially. Bifurcations from TW are Hopf ones, and for WR it is pitchfork. 

The spatial structure of magnetic fields is shown on Fig. [11] A common dominant feature in the majority 
of nonlinear convective hydromagnetic regimes is concentration of magnetic field flux in half-ropes located 
near the horizontal boundaries. In the primary attractors this structure is inherited from the respective 
kinematic dynamo modes. Concentration of magnetic field near boundaries of the layer was observed by St 
Pierre (1993) in simulations of subcritical magnetic field generation by thermal convection of rapidly rotating 
fluid (for the same boundary conditions, as employed here). Such behaviour of magnetic field was observed 
in many computations (see a discussion and references in Zheligovsky, 2010) and it is usually expected for 
electrically perfectly conducting boundaries (although Zheligovsky, 2010, found for these boundary condi- 
tions an example of a time-periodic nonlinear convective magnetic dynamo, in which magnetic field always 
remains concentrated inside the layer of fluid). 

In what follows we overview the MHD attractors found in computations and their bifurcations. 
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Figure 11: Isosurfaccs of magnetic energy density of magnetic fields, at the level of a half of the maximum, in steady convective 
MHD attractors: Ta = 150, Sf ^ (a); Ta = 200, S^^ (b); Ta = 250, Sj^ (c); Ta = 300, Sj^ (d); Ta = 430, Sf^ (e); Ta = 600, 
Sf^ (f); Ta = 675, S^^ (g); Ta = 680, Sf^ (h); Ta = 721, S^^^- (i);'Ta = 690, Sf° (k); Ta = 740, Sf° (ij; Ta = 775, Sf° 
(m); Ta = 1175, Sj^ (n). Stagnation points of the flow on the horizontal boundaries are shown by dots. Four periodicity cells 
are displayed. 
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Figure 12: Bifurcation diagram of the MHD system for < Ta < 200. Solid lines denote stable branches, dashed unstable, 
thin dashed conjectured. 



6.1. MHD attractors emerging from TW, < Ta < 81 

The bifurcation diagram of MHD regimes at the interval < Ta < 200 is shown on Fig. [T21 
The interval of Ta, where TW exists, consists of two subintervals of kinematic dynamo action, separated 
by a window 82 < Ta < 86 of non-generating TW regimes. Growing magnetic modes generated in the 
two subintervals differ, for instance, in their symmetries (see Table [2]). This results in emergence of MHD 
attractors of different types. At the lower subinterval of Ta the primary MHD attractor is Pj''^, similarly 
to TW drifting along both horizontal axes, x and y, and periodic in the co-moving reference frame. At the 
right end the MHD branch Pj""^ terminates on TW at Ta — 81.80, where the respective magnetic eigenvalue 
of TW crosses the imaginary axis; to the left it continues to Ta = 0. 

6.2. MHD attractors emerging from TW, mode interaction 

As noted in Section 3, TW bifurcates supercritically in a Hopf bifurcation from Ri as Ta decreases below 
188. For the employed value of Pm, at the point of bifurcation Ri generates magnetic field kinematically. 
The respective magnetic eigenvalue is real and not large, hence attractors observed for Ta not far from 188 
can be related to interaction of two instability modes of Ri, hydrodynamic and magnetic ones. 

The symmetry group of Ri is generated by S2, ''72/2' 7^ ^^'^ 9- hydrodynamic eigenspace the 

symmetries S2 and 7^ act non-trivially, and in the magnetic eigenspace the action of the symmetries S2 , q 
and 7^ is non-trivial. Hence, the results of Golubitsky et al. (1988) (Section XX §2) on Hopf / steady-state 
mode interaction with the 0(2) symmetry group are applicable. (The symmetry q plays no role, since, in 
the notation ibid., q = (7r,7r) is a superposition of two shifts, in = tt and = 7:, which are elements 
of SO(2) and S^, respectively. In our problem, the action of SO(2) is generated by the shifts 7^ along 
the direction of the rolls.) According to Golubitsky et al. (1988), the trivial steady state (Ri in our case) 
bifurcates with emergence of an S2-symmetric steady state, when a real eigenvalue becomes positive, and 
with a simultaneous emergence of rotating (i.e. travelling, in our parlance), TW, and standing waves, SW, 
when a complex pair of eigenvalues crosses the imaginary axis. The rotating wave can further bifurcate to 
a modulated rotating wave (which in fact is a 2-torus). 

In our system, in agreement with this general theory, Ri bifurcates to TW (in the hydrodynamic sub- 
space) and to 82^^. TW further bifurcates to a modulated travelling wave, Pj^ (which we classify as a 
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Tabic 5; Families of attractors F^* and C^^. Column 2 presents the type of attractors (disregarding drift frequencies); 
P periodic, Q quasiperiodic, C chaotic, column 3 basic frequencies of periodic and quasiperiodic regimes, columns 5 and 6 
time-averaged kinetic and magnetic energies, respectively, the last column the list of individual runs as pairs ra(integration 
time) . 
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periodic orbit, since we omit drift frequencies in the description of attractors). At Ta = 87 the bifurcation 
where Pj^ emerges is subcritical, and close to the point of bifurcation the periodic orbit is unstable. It 
gains stability as it turns back in a saddle- node bifurcation at Ta = 36. We classify the stable periodic 
orbit observed in computations as the modulated rotating wave predicted by the theory, because 1) of the 
similarity of the spatial structure of TW and the flow in Pj^, as well as of the dominant magnetic mode 
of TW and magnetic field in Pj^; 2) the symmetry group of Pj^ is the one expected on the theoretical 
grounds for the branch emerging at Ta — 87 from TW due to the mode interaction (note that it includes a 
spatio-temporal symmetry: the symmetry about a vertical axis with a shift in time by a half of a period); 
and 3) the temporal frequency of Pj^ is close to the Hopf frequency of TW. 

The steady state Sf-^ bifurcates as Ta is decreased to a less symmetric steady state Sf ^ , which is outside 
the Hopf / steady-state mode interaction center manifold. As Ta is decreased further, S^^ becomes unstable 
in a Hopf bifurcation with a stable periodic orbit emerging. We label it P^^, because its spatial structure and 
temporal frequency are similar to those of the unstable standing wave SW emerging from Ri simultaneously 
with the travelling wave TW, and its symmetry group is a subgroup of the one of SW. The general theory of 
Hopf / steady-state mode interaction predicts two types of periodic orbits bifurcating from SW, but judging 
by its symmetries none of them is our P^^. The conjectured relation of P^w ^fj]^]^ jg shown by thin 
dashed lines on Fig. [T^ 

6.3. MHD attractors emerging from TW, 36 < Ta < 127 

Both periodic orbits, Pj^ and P^w^ gjyg j.jgg complex families of attractors, F"""^ and C^^, discussed 
in this subsection (see Table [5]). Like Pj^, all attractors constituting F"^^ drift along both horizontal axes 
and thus have two drift frequencies. The family F"""^ starts from the primary periodic MHD attractor Pj^ 
at Ta = 36. 

As Ta is increased beyond Ta w 60, the second frequency, /2, appears in a Hopf bifurcation, i.e. a stable 
torus emerges. Afterwards, a sequence of bifurcations of halving of /2 takes place (often such a sequence 
is called a cascade of period doubling bifurcations for a torus; we prefer to abstain from the use of this 
terminology, since the concept of a period applied to a torus is not all too transparent). This sequence of 
bifurcations is analogous to the Feigenbaum (1978) scenario of period doublings for a logistic map, and we 
label this family of attractors by F. 

Logistic map is defined by the recurrent relation Xn+i = rxni^ ~ ^n), < x < 1. When r is increased 
over 3, a stable 2-cycle is created in a flip bifurcation. Subsequent flip bifurcations result in emergence of 
cycles of lengths 4, 8, . . .. The bifurcations accumulate at r = 3.57. For larger r, chaotic behaviour sets in, 
alternating with windows of periodic orbits of periods (2m -I- 1)2*^. For r > 4, x„ escapes from the interval 
< X < 1 and diverges. 

In the convective hydromagnetic system, on increasing Ta we have detected three halvings of the fre- 
quency /2 ~ 0.15, followed by a chaotic behaviour. For larger Ta, lower basic frequencies set in, /2/3 for 
75 < Ta < 77 and /2/6 for 78 < Ta < 78.5; their appearance agrees with the general theory. No more 
halvings have been observed. This is consistent with the fact that a complete sequence of period doubling 
bifurcations of tori is a structurally unstable scenario (see Coullet, 1984), and generically it is interrupted 
by the onset of chaos. (We could also miss some halvings just having considered not enough values of Ta.) 
For still larger Ta, a frequency /2 ~ /2/2 is observed. On the one hand, emergence of this quasiperiodic 
regime is clearly outside the framework of the Feigenbaum period doubling scenario (cf. Fig. [13] (b) and 
(h) ), suggesting that attractors for 79 < Ta < 80 do not belong to our Feigenbaum family. On the other, 
the standard theory takes into account only two first terms in the Taylor expansion of the dynamical system 
on the central manifold, and terms of higher order may be responsible for the birth of this new regime. 

Plots of the energies (Fig. [T5)) . not affected by drift frequencies, clearly display the frequency halving. 
By contrast, plots of the time dependencies of real or imaginary parts of individual Fourier coefficients (not 
shown) are marred by the extra drift frequencies, making these plots by far more obscure. A similar effect 
is observed in Poincare sections. Frequency halving for tori usually shows itself unambiguously on Poincare 
sections, however, the standard Poincare sections are not particularly enlightening when additional drift 
frequencies are present. Their infiucnce is eliminated in Poincare sections, which are constructed with the 
use of absolute values of Fourier coefficients (cf. Figs. [14] and [15]). 
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(g) (h) 

Figure 14: Poincarc sections of regimes in the family F''"^ on the (Rcfij^j 2 ijRcS^j j 1) plane (horizontal and vertical axes, 
respectively) defined by the condition Re{)li_2,i =0- ^» = 65 (a), Ta = 70 (b), Ta = 71.5 (c), Ta = 72 (d), Ta = 74 (e), 
Ta = 77 (f), Ta = 78.5 (g), Ta = 80 (h). 
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The family C exists for 78 < Ta < 126. It is comprised of chaotic attractors, alternating with windows 
of periodic and quasiperiodic regimes (Fig. I16p . The family starts with the periodic regime psw^ which is 
stable for 120 < Ta < 126. The orbit possesses a symmetry, which is a combination of S2 with a shift by 
a half of its temporal period, hence, in agreement with Krupa (1990), the periodic orbit does not drift. As 
Ta is decreased, the symmetry is lost, the emerging non-symmetric orbit and subsequent attractors of the 
family have two drifting frequencies (not discussed here). Subsequently, a second frequency appears in a 
Hopf bifurcation, which is halved afterwards. On a further decrease of Ta, we observe an intermittency 
of chaotic and quasiperiodic attractors. The sequence "a periodic orbit, a torus, chaos" is standard; it is 
usually explained by appearance of the third basic frequency, which makes the system structurally unstable 
and results in the onset of a chaotic behaviour (Ruelle and Takens, 1971). The observed windows of 
quasiperiodicity can be attributed to frequency locking (see, e.g., Ott, 2002). 

Notably, coexistence of attractors of different types is observed: two types of MHD attractors (F"^^ and 
Pj^) in the interval 36 < Ta < 81, and a (magnetically stable) hydrodynamic attractor (TW) with an MHD 
attractor (C^^) in the interval 81.80 < Ta < 87. Coexistence of a hydrodynamic and MHD attractors can 
be described in physical terms, as stiff excitation of a magnetic field: small magnetic perturbations of TW 
decay to the hydrodynamic state, while large magnetic perturbations give rise to MHD regimes. Furthermore, 
in the interval 78 < Ta < 80 three MHD attractors coexist: Pj"^ , F^"" and C^^. From the point of view 
of the theory of dynamical systems there is nothing extraordinary in coexistence of attractors in a nonlinear 
system, however, to the best of our knowledge coexistence of three convective MHD attractors was never 
observed before. 
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Figure 17: Phase portrait of the MHD system on the interval 216 < Ta < 224. 



6.4- MHD attractors emerging from Ri, 127 < Ta < 506 

Dominant magnetic modes of Ri are of the same type for all Ta. As attested by symmetries, the primary 
MHD attractor is the steady state S^^ detaching from Ri at Ta = 725.3 in a subcritical pitchfork bifurcation. 
It exists (i.e. the MHD states are stable) for 170 < Ta < 216. On decreasing Ta, a steady state S^^ with 
a smaller group of symmetries emanates in a pitchfork bifurcation; it becomes unstable and bifurcates to 
in a Hopf bifurcation. 

On increasing Ta, in a small interval 216 < Ta < 224 four bifurcations occurring in the vicinity of 82'^ (in 
the phase space) are observed (see Fig.[7](b) ). They take place in two distinct invariant subspaces, Fix(rs2) 
and Fix(s2) (Fix(s) denotes the set of fixed points under the action of a symmetry s). Developments in each 
invariant subspace are mutually independent (see a sketch of geometry of the phase space shown on Fig. 1171 
the vertical plane represents Fix(rs2) and the horizontal one Fix(s2) )• 

The first (on increasing Ta) one is a Hopf bifurcation in Fix(rs2) at Ta = 216.13, in which a periodic 
orbit, P^^^, emanates. Another one is a saddle-node bifurcation at Ta = 223.1 in the same invariant subspace, 
in which a branch of steady states, 84^, terminates. Near the point of bifurcation the basin of attraction 
of the steady states becomes vanishingly small, requiring short steps when continuing the branch in Ta 
so that computed trajectories were not attracted by P^^^. On decreasing Ta from, say, 224, eigenvalues of 
stability modes are in the beginning imaginary, but both the real and imaginary parts decrease in magnitude, 
apparently tending to zero. At Ta ~ Tac, 223.101 < Tac < 223.10075, the discriminant of the quadratic 
characteristic equation, defining the two eigenvalues, changes its sign. As a result, the dominant eigenvalue 
exhibits a counterintuitive behaviour: Coefficients of the characteristic equation depend almost linearly on 
Ta near this value Tac, but the change of sign of the discriminant implies, that the dominant real eigenvalue 
is continuous, but non-smooth (see Fig. \TE\\ . The graph of the growth rate of the dominant mode bends at 
Tac and for Ta < Tac behaves like \/Tac — Ta. The largest real eigenvalue starts to grow much faster, and 
a saddle-node bifurcation occurs between Ta = 223.09999 and 223.1. An almost simultaneous vanishing 
of both real and imaginary parts of eigenvalues implies, that two parameters are necessary to describe this 
bifurcation, i.e. we are in a vicinity of a codimension two bifurcation. Variation of Ta near the left end of 
the branch S^^ is equivalent to a motion along an one-dimensional curve on the plane of the parameters. 

Vanishing of a pair of complex eigenvalues of linearisation of a dynamical system is an attribute of the 
Takens-Bogdanov bifurcation (see Guckenheimer and Holmes, 1988). There are further indications that 
the bifurcation occurring in Fix(rs2) close (in the parameter space) to the values, which are fixed in our 
simulations, might be a Takens-Bogdanov bifurcation: in this interval, 216 < Ta < 224, a periodic orbit 
emerges and afterwards becomes unstable or disappears and a steady state emerges, and they are close 
neighbours in the same symmetric subspace. However, a more attentive inspection of the system convinces 
that this conjecture is wrong. 

Only ^ might be the trivial steady state suffering the bifurcation (although the behaviour of the 
eigenvalues reminiscent of the Takens-Bogdanov bifurcation is registered for S , this branch turns back in 
a saddle-node bifurcation and thus can not serve as the trivial steady state, existing in the case of Takens- 
Bogdanov bifurcation for all parameter values in the vicinity of the point of bifurcation). It is quite possible 
that a pair of complex eigenvalues of linearisation of 82^ simultaneously become zero upon a variation 
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Figure 18: Computed mean of the two dominant eigenvalues (blaek circles; coincides with the real part of the eigenvalues for 
Ta > Tac ~ 223.101, where the two eigenvalues are complex conjugate) and the discriminant of the quadratic characteristic 
equation (open circles) for the dominant eigenvalues of the operator of linearisation versus Ta (horizontal axis) near the left 
end of the interval of existence of 84^^^. Fade extrapolation of the mean of the two eigenvalues (solid line), the discriminant of 
the characteristic equation (dotted line) and the resultant extrapolated dominant real eigenvalue for Ta < Tac (dashed line). 

of Ta and a second parameter: 82^^ does have small in magnitude complex eigenvalues in the interval 
216 < Ta < 220 for the parameter values that we have employed; also it is not far from Sj^, suggestive of a 
similar behaviour of eigenvalues of the two branches. The structure of our system is compatible only with the 
Takens-Bogdanov bifurcation (with the Z2 symmetry group generated by S2, which S^^ does not possess) 
involving a stable periodic orbit and stable steady states, distinct from the trivial one; such a diagram is 
shown on Figs. 7.3.7 and 7.3.9 ibid. On these Figures, the trivial steady state (the analogue of 82^) is stable 
for /ii < 0, /i2 < (using the parameter notation ibid.), the other two steady states (the analogues of 84^ 
and its symmetric counterpart) for fii > max(0,/i2), and a periodic orbit for fi2 > max(0,c/ii), c ~ 0.752. 
Hence variation of Ta on the interval 216 < Ta < 224 must be equivalent to a motion along a curve on the 
(/^ii/^2) plane, beginning in the quadrant /ii < 0, /i2 < 0, following to the quadrant /ii < 0, /i2 > 0, and 
passing to the region /ii > 0, fJ.2 < c/ii. Consequently, the curve must cross the line fj,i = 0, where in the 
Takens-Bogdanov bifurcation a branch of steady states (analogues of 84^), stable or unstable depending on 
the sign of /ii, emanates from the trivial steady state in a pitchfork bifurcation. This is inconsistent with 
the behaviour of eigenvalues of linearisation of 83^ (we have computed these unstable steady states after 
the Hopf bifurcation, giving rise to P^^, by imposing the symmetries rj^^^ and S2) - the dominant and 
subdominant eigenvalues do become real between Ta — 221 and 221.5, but they remain positive and do not 
become small on the interval 221.5 < Ta < 224. Therefore we conclude, that bifurcations in our system are 
not induced by a Takens-Bogdanov bifurcation (a finite perturbation of a bifurcation of this type, in which 
more than two parameters are essential for the description of the bifurcation, is not ruled out). Further 
analysis is necessary to understand it theoretically. 

At Ta — 216.75 an eigenvalue of linearisation of (1) near 8^^ with the associated eigenvector in Fix(s2) 
(the dominant eigenvalue of the restriction of the operator of linearisation to this subspace) becomes positive. 
In this supercritical pitchfork bifurcation a steady state 8^^"'^ emerges. 8ince at the bifurcation 8^^ is unstable 
with respect to perturbations in Fix(rs2), the branching steady state is also unstable to such perturbations 
near the point of bifurcation. 8^^ gains stability at Ta — 218.64 in a subcritical Hopf bifurcation. The 
branch is stable up to a subcritical pitchfork bifurcation at Ta = 271.66. 8^^ is the only MHD steady state 
that we have found, whose group of symmetries does not involve superpositions of the symmetry r with any 
other symmetries. As a result, magnetic patterns on top and bottom of the layer are different (cf. Fig. [TT] 
(c) and other panels on this figure). 

We again encounter coexistence of three distinct MHD attractors in the interval 223.11 < Ta < 223.3: 
83^, Sf^ and P^^. The steady state 8^-^ becomes unstable at Ta — 466.5 in a subcritical pitchfork bifurca- 
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Figure 19; Phase portrait of the MHD system near Ta = Ta^ ~ 506.07 . 

tion. 

The symmetry group of 85^ is a subgroup of the one of 82^^ (see Table [3]). That the steady states are 
related and some symmetries are lacking, is seen on Fig. [11] (b) and (e) displaying magnetic patterns of the 
steady states. These facts suggest that S^"'^ bifurcates from the unstable Sf-^. At both ends of the interval of 
stability, S^^ undergoes supercritical Hopf bifurcations with stable periodic orbits, Pf-^ and Pf^^, emerging. 
The periodic orbit P^^ with the same group of symmetries as that of 85^ is observed for 343 < Ta < 377. 

The other periodic orbit bifurcating from S^^, Pf'^, exists in a small interval 505.1 < Ta < 506 and it 
has a much smaller symmetry group than 85^^ . 

6. 5. MHD attractors emerging from Ri , Ta near 506 

At the right end of the interval of existence of the periodic orbit P3 ^ , it terminates near Ta — 506 on a 
steady state 8g^"'^ in a bifurcation, which appears to be similar to a saddle-node bifurcation on invariant circle 
(8NIC; see Izhikevich, 2006). However, since the symmetry group of the system is non-trivial, the details 
of the bifurcation in our case differ from those of the canonical 8NIC. In a generic system, a periodic orbit 
emanates in a 8NIC bifurcation subsequent upon a saddle-node bifurcation of a steady state. The 8NIC 
occurs under the condition that for any parameter value before the saddle-node bifurcation (we only consider 
a small neighbourhood of the critical value of the bifurcation parameter), both parts of the one-dimensional 
unstable manifold of the unstable steady state terminate on the stable one. At the point of bifurcation, 
the two steady states collide and a homoclinic trajectory emerges. After the bifurcation, a periodic orbit is 
formed from this homoclinic trajectory. 

8NIC is generic in one-parameter dynamical systems and it is often observed in simulations, where a 
parameter is varied. At Ta = Tac ~ 506.07 the periodic orbit Pg^ terminates on a heteroclinic cycle; 
occurrence of this upon variation of just one scalar parameter requires some degeneracy of the system, since 
a heteroclinic (homoclinic) connection from an equilibrium to an unstable equilibrium in a generic system 
occurs only for a singular parameter value (such connections are structurally unstable) . For formation of a 
heteroclinic cycle, several such connections must happen simultaneously. 

Existence of a heteroclinic orbit at Ta = Tac becomes possible due to the presence of a non-trivial 
symmetry group. 8tructural stability of homoclinic and heteroclinic connections to an unstable equilibrium 
in symmetric systems relies on the presence of symmetry-invariant subspaces (Guckenheimer and Holmes, 
1988). A sketch of geometry of the phase space of our system is shown on Fig. [191 The steady state 8^^ 
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Figure 20: Kinetic (upper curve, left vertical axis) and magnetic (lower curve, right vertical axis) energies for Ta = 506, P^^ 
as a function of time (horizontal axis) for a trajectory starting at the steady state S^^, the MHD attractor for Ta = 505. 

possesses a symmetry group isomorphic to D4 k Z2. The plane represents the fixed point subspace for 
the steady state symmetry group, the vertical direction the antisymmetric complement. We describe the 
bifurcation starting from larger values of Ta slightly above the saddle-node bifurcation of S^^. Denote by 
the unstable counterpart of Sg^; it belongs to Fix(D4 k Z2). The one-dimensional unstable manifold of 
S^^ also belongs to this subspace. A part of the unstable manifold terminates on the stable Sg^, another 
one on Ri (this has been checked numerically), which is stable within this subspace. The connections from 
S^*'"'^ to Sg^ and Ri are structurally stable. The one-dimensional unstable manifold of Ri terminates on the 
stable S^^ (we have also checked this numerically). 

At Ta = Tac Sf ^ and Sg^^ collide, thus creating the structurally unstable heteroclinic cycle Sg^ — Ri — ?> 
Sg'"'^ (it exists only for a single value of Ta). For Ta < Tuc, S^^ disappears and we observe a periodic orbit 
in place of the heteroclinic cycle. The period of this orbit tends to infinity, as Ta approaches the point of 
bifurcation Tac. 

Kinetic and magnetic energies for the periodic orbit for Ta = 506, which is close to the point of the 
saddle-node bifurcation (hence the period of the orbit is large) is shown on Fig. 1201 The minima of magnetic 
energy are attained when the trajectory is close to Ri; from these steady states it jumps rapidly to the 
former Sg^ (kinetic energy becoming close to 161), from where it slowly moves back towards Ri. 

Notably, close to the critical Ta we observe not a periodic, but rather a chaotic behaviour (see Fig. [20)1 . 
The non-periodicity (the behaviour is, loosely speaking, periodic, but "periods" significantly vary) can be 
caused by reasons of numerical nature, such as round-off errors or other numerical noise accumulating during 
the long periods. As noticed by Busse and Heikes (1980), round-off errors are comparable with the amplitude 
of the integrated field near stagnation point, significantly affecting trajectories approaching a heteroclinic 
cycle: Numerical noise can result in emergence of quasi-periodic regimes with randomly varying "periods" 
(Stone and Holmes, 1990), like we observe here, as well as strictly periodic cycles of a large period (Nore et 
al, 2003), or alter the amplitude of cycles (Stone and Armbruster, 1999). Numerical noise can be expected 
to affect similarly large-period orbits near the critical parameter value. 

6.6. MHD attractors emerging from Ri, 506 < Ta < 718 

The primary steady state Sg^^, connected with the S^^ branch by the periodic orbit Pf^^, is different 
from Sf-^ and other S^^ steady states considered above in that it is related not to the dominant but a 
subdominant magnetic mode of Ri . The symmetry group of the steady state S^^ is not a subgroup of the 
symmetry group of S^^ (coinciding with the symmetry group of the dominant magnetic mode) and it has a 
twice smaller period in the direction along the axis of rolls, indicating that S^^ is unrelated to the dominant 
magnetic mode (cf. Fig. [S] (c) and Fig. [TT] (f) ). We have computed the subdominant magnetic mode by 
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restricting the kinematic dynamo problem for Ri on the subspace Fix(g7^^^). The subdominant magnetic 

mode, which is dominant in Fix(g7^^^), has the same group of symmetries as that of S^^ and a spatial 

structure similar to that of the magnetic field of Sg^^ (see Fig. [11] (f) ). The associated eigenvalue is real and 
positive in the interval 509 < Ta < 682 and admits the maximum 0.4 at Ta — 598. On increasing Ta, Sg^ 
becomes unstable in a pitchfork bifurcation with emergence of the steady state Sy^, becoming unstable in 
its turn in a subcritical pitchfork bifurcation at Ta — 682.26. 

The symmetry group of S^*^^ is a subgroup of the symmetry group of S^^, indicating that the former 
state has possibly bifurcated from the latter one. On decreasing Ta, Sf^ becomes unstable in a subcritical 
pitchfork bifurcation; in the interval 674 < Ta < 682 the steady state coexists with Sy^. When Ta is 
increased, a stable periodic orbit emanates in a Hopf bifurcation from Sg^. On increasing Ta further, 
this periodic orbit terminates on a structurally unstable heteroclinic cycle. We discuss this bifurcation in 
detail in the following subsection. 

6.7. MHD attractors emerging from Ri, Ta near 718 

We have discussed in subsection 4.3 a periodic orbit terminating on a structurally unstable heteroclinic 
cycle involving connections between two steady states. At Ta = Tac ~ 718.16 a similar bifurcation takes 
place. The periodic orbit Pj^ terminates on a structurally unstable heteroclinic cycle with three steady 
states involved, when one of the steady states undergoes a saddle-node bifurcation. 
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A sketch of geometry of the phase space is shown on Fig. 1211 The plane represents the hydrodynamic 
subspace invariant under the symmetry q, and the vertical direction the complementary magnetic subspace. 
We describe the bifurcation starting from the larger values of Ta slightly above the point Ta = Tqc of 
the saddle-node bifurcation of WR. For such Ta, in the hydrodynamic subspace we observe Ri, stable in 
the hydrodynamic subspace, stable WR and its unstable counterpart, WR. A part of the one-dimensional 
unstable manifold of WR terminates on Ri, the other one on WR (since WR emerges in a subcritical 
bifurcation from Ri and it disappears in a saddle- node collision with WR) . Connections from WR to Ri and 
WR are structurally stable, because Ri and WR are stable in the hydrodynamic subspace. The connection 
from Ri to belongs to Fiyi{q'-f^^^) and it is structurally stable, because S^^ is stable in this subspace 
(this has been checked numerically). The unstable manifold of 82^^ belongs to Fix(s2) and it terminates on 
WR, stable in this subspace. 

At Ta = Tac, WR and WR collide creating the heteroclinic cycle WR Ri ^ S^^ WR in the 
subspace Fix(s2)- The cycle is asymptotically stable within this subspace, but not in the entire phase 
space, because WR possesses a growing magnetic mode. The cycle is structurally unstable; it exists only 
for Ta = Tac- For Ta < Tac, WR disappears and we observe a periodic orbit in place of the heteroclinic 
cycle. The period of the orbit tends to infinity, as Ta approaches the point of bifurcation. Surprisingly, the 
orbit is asymptotically stable, despite it has bifurcated from an asymptotically unstable heteroclinic cycle. 

We illustrate the behaviour described above by plots of kinetic and magnetic energies and discrepancies 
for the symmetries S2 and 97^/2 ^ periodic orbit at Ta = 717.9 (Fig. [22t . near the point of bifurcation, 
where the orbit is similar to the structurally unstable heteroclinic cycle. Plateaux of constant values of 
kinetic energy on Fig. [22] (such as 1870 < t < 1930 and 1990 < t < 2060) represent time intervals when 
the trajectory is near the steady states Ri and S^^ , respectively. The inflection point near t = 2090 shows 
where the trajectory is in the vicinity of the former steady state WR. Magnetic energy is small during the 
transition from the former WR to Ri, attaining a minimum at t — 2110. The symmetry present 

during the transition from Ri to S^^. An exponential growth of the symmetry discrepancy takes place 
during the departure from Sj^. The symmetry S2 is slightly broken near the former WR. 

The discrepancy of the symmetry S2 is shown for an axis changing position in time (the position is 
optimised to minimise the discrepancy). The y coordinate of the moving axis is plotted on Fig. [25] (c); the 
displacement in the x direction is by several orders of magnitude smaller, and it is described by a function 
of a similar shape. The motion of the axis can be regarded as a drift of the cycle Pf^^ along a group orbit. In 
accordance with the theory of Krupa (1990), the axis is steady when the symmetry S2 is present; it moves, 
when the discrepancy is maximal (for instance, at the time interval 2140 < t < 2160). 

For Ta = 719 a trajectory starting near S^^ visits the following steady states: 82^ WR Rd —>■ 

6.8. MHD attractors emerging from WR, 721 < Ta < 725 

The primary steady state emerges from WR in a pitchfork bifurcation at Ta ss 725.71, when the 
magnetic growth rate becomes positive on decreasing Ta. It is stable in a short interval of Ta and disappears 
in a saddle-node bifurcation at Ta « 720.25. 

6.9. MHD attractors emerging from Rd, 685 < Ta < 787 

We have found numerically that Rd possesses two types of dominant magnetic modes (see Table [1] and 
Fig. [3]). In the interval of Ta where Rd exists there are two windows of kinematic dynamo action, resulting 
in two intervals of nonlinear dynamos; the lower one is 684 < Ta < 787. 

The first (on increasing Ta) attractor S^^^^ is primary, but it is related to a subdominant magnetic mode. 
The mode has the same spatial structure as the respective nonlinear steady state displayed on Fig. [TT] (k). 
The mode has the period L/3 along the axes of rolls (in contrast with periods L/2 and L for two other 
magnetic modes. Fig. [TT] (m) and (n) ). The associated eigenvalue is real and positive for 685 < Ta < 711, 
the maximum is 0.1 . 
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Figure 22: Kinetic (solid line) and magnetic (dashed line) energies (a), symmetry discrepancies (b) and tiie y coordinate 
of the axis of the symmetry S2 (c) versus time (liorizontal axis) for Ta = 717.9, Pj'^. On (b), tiic following symmetry 
discrepancies are traced: in the flow, 7]^y2 (solid line), i'l^/2 (dotted line), S2 (dash-and-three-dots line); in magnetic fleld: 
'^L/2l (short-dash line), rjf^^^ (dash-and-dot line), S2 (long-dash line). Discrepancy of a symmetry s in a field f is measured 

as f \{ — s(f)P<ix/ / Ifpdx, where integration over a periodicity cell is assumed. 
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Another primary MHD attractor is bifurcating from RD in a pitchfork bifurcation when the respec- 
tive eigenvalue of the kinematic dynamo problem crosses the imaginary axis. The two branches of primary 
steady states are connected by the branch S^'^, with a smaller symmetry group (cf. Fig. [Tl](k), (1), (m) ). 

6.10. MHD attractors emerging from Rd, 1118 < Ta < 1355 

The third primary branch emanating from Rd, the steady state 84^^, is always stable when exists. Both 
ends of the branch terminate on Rd at Ta = 1117.47 and Ta = 1355.20 . 

7. Conclusion 

Our results demonstrate that the influence of rotation on magnetic field generation by thermal convection 
is non-monotonic and in no way simple. Its nature can only be understood by a careful investigation of 
attractors of the underlying dynamical system and bifurcations of branches of convective MHD regimes. 
Such an investigation is made more exciting by the presence of a large group of symmetries of the dynamical 
system. Even close to the critical point for the onset of convective motion, in the terms of the Taylor 
number, the bifurcation diagram is quite complex. There are several intervals of coexistence of two or three 
convective MHD attractors (also sometimes with a non-generating hydrodynamic one). 

An overall picture that we have obtained fits well the usual beliefs, that up to a point an increase of 
the rate of rotation from zero benefits magnetic field generation both in linear and nonlinear regimes, and 
after attaining a maximum of the mean magnetic energy (in our simulations, Em = 32.30 in the regime S^^ 
for Ta — 216) on a further increase of the rotation rate _E,„ gradually falls off till magnetic field generation 
ceases, and at still higher angular velocities the fluid flow is arrested. However, the presence of many 
individual branches and windows of nonlinear dynamo action (e.g., between the MHD steady states 83^^ 
and 84^^) adds a rich "small-scale structure" to this otherwise "smooth" general picture, on some intervals 
of Ta even reversing it. For instance, an increase of the rate of rotation from zero inhibits magnetic field 
generation in the periodic regime P"^^, rather than enhances it. 

Apparently the conjecture, that increasing the Taylor number first helps, afterwards hinders and finally 
halts magnetic field generation, has not yet been proved, but it appears to be correct for all parameter values 
(provided the Rayleigh number is sufficiently large). In the absence of magnetic field the critical Rayleigh 
number for the onset of convection grows with the Taylor number (see Chandrasekhar, 1961); thus, in the 
hydrodynamic case fast rotation is guaranteed to arrest convective flows. If a magnetic field is present, the 
system is more complex and various phenomena — such as subcritical loss of stability of the basic state with 
the fluid at rest — are not ruled out for sufficiently strong initial magnetic fields; this can result in extension 
of the interval of the Taylor number, for which convective MHD states persist for the given other parameter 
values. However, it is unlikely that this interval becomes infinite even for Pm S> P, since — in physical 
terms — magnetic field generation relies on advection by the fiuid, and not vice versa. This question awaits 
further investigation. 

We have observed a number of interesting bifurcations. To the best of our knowledge, two global 
bifurcations were not observed before. They are similar to the 8NIC (saddle-node on invariant circle) 
bifurcation, and their peculiarity stems from the presence of a non-trivial symmetry group in the convective 
MHD system. In these bifurcations, at Ta w 506.07 and Ta 718.16, a periodic orbit terminates not 
on a homoclinic (as in the 8NIC), but a heteroclinic cycle, whose existence relies on the presence of the 
symmetry group. Among more common, albeit rather seldom observed in natural systems families, that 
we have encountered, an incomplete Feigenbaum sequence of "period doubling bifurcations of a torus" , 
pTW (occurring between Ta = 57 and 80) is notable, as well as an intermittent sequence of chaotic and 
quasiperiodic regimes, (occurring for 78 < Ta < 126), in which existence of the quasiperiodic regimes 
can be apparently linked with frequency locking. 

Another unusual set of bifurcations observed on the interval 216 < Ta < 224 is associated with a pair 
of complex eigenvalues of the linearisation of the dynamical system, which are simultaneously vanishing. 
Numerical examination of eigenvalues of the supposed trivial steady state 82^ suggests that the observed 
bifurcations in our convective MHD system are not linked with a Takens-Bogdanov bifurcation occurring 
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for nearby parameter values, but we do not rule out a finite perturbation of a bifurcation of this type. We 
plan to perform a further investigation of the bifurcation in order to gain a mathematical understanding of 
its nature. 

We have found a number of branches of MHD steady states, which are parity-invariant (i.e. have the 
symmetry rs2) or possessing the symmetry about a vertical axis, S2. Certain periodic orbits (Ta = 57, the 
family F™; Ta = 123, the family 0*=*^; 216.13 < Ta < 223.3, branch Pf^) are not pointwise symmetric, 
but symmetric on average, i.e. the symmetry with a time shift by a half of temporal period is present. 
Possession of such a symmetry is an important property of attractors, since in convective MHD systems 
with this symmetry the global a— effect is zero. Whilst in the presence of the global a— effect the system 
is inherently unstable to large-scale perturbations, in its absence the instability, when present, develops on 
time scales of a larger order. In the latter case it is described by a highly complex nonlinear mixed system 
of PDE's of the second and third order (see Zhcligovsky, 2009b), incorporating such physical effects, as 
combined eddy diffusivity and eddy advection. We are planning to follow this line of research. 

For the chosen values of the Rayleigh and magnetic Prandtl numbers, which are not particularly high, the 
convective system is not far from the onset of convection and magnetic field generation, and the collection of 
hydrodynamic and magnetic structures that we encounter is not rich (in the MHD regimes the flows take the 
form of perturbed rolls, and magnetic field concentrates in half-ropes located near the boundaries). More 
complex structures are produced in a more vigorous convection for higher values of the Rayleigh and the 
magnetic Prandtl numbers. 

A natural extension of our investigation consists of continuation of branches of regimes, that we have 
found, in other parameters, which would significantly contribute for gaining an insight into the mathematics 
and physics of convective dynamos (i.e., of the geometry of branches of attractors in the parameter space). 
This project is, obviously, quite resource demanding, but in our opinion this approach is more promising 
for the deeper understanding of convective MHD systems, than analysis of individual runs for extreme 
parameter values. 
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